C""-THEORY OF DAMPED WAVE EQUATIONS WITH STABILISATION 



FUMIHIKO HIROSAWA AND JENS WIRTH 



Abstract. The aim of this note is to extend the energy decay estimates from [J. Wirth, 
Wave equations with time- dependent- dissipation I. Non-effective dissipation., J. Differential 
Equations 222 (2006) 487—514] to a broader class of time-dependent dissipation including very 
fast oscillations. This is achieved using stabilisation conditions on the coefficient in the spirit of 
[F. Hirosawa, On the asymptotic behavior of the energy for wave equations with time- depending 
coefficients, Math. Ann. 339/4 (2007) 819-839]. 



1. The problem under investigation 

We want to investigate the Cauchy problem 

□u + 26(i)ut = 0, u(0, •) = ui, Dtu(0,-) = U2 

for a weakly damped wave equation with time-dependent dissipation, as usual □ = 
denotes the d'Alembertian and D = ^id. For this we apply a partial Fourier transform to get 
the ordinary differential equation 

utt + \^\'^u + 2b{t)ut^0 

parameterised by the frequency ^. To formulate a first order system corresponding to this second 
order equation we consider V — Dtu)'^, such that 

^^'^^(lei 2il%)^-^(*'^)^- 
We denote its fundamental by £{t, s,^), i.e. 

Dt£{t, s, = A{t, 0£{t, s, 0, - ^ e 

Our aim is to understand the structure of £(t, s,^) and its asymptotic behaviour as t — s- oo in 
dependence on the coefficient function b — b{t) extending results from [9J. If b{t) = vanishes 
identically, we denote the fundamental solution as £o(t,s,S,) and refer to it as free solution. 

There is a very strong interrelation between properties of the coefficient function b — b{t) and de- 
cay properties of solutions to the above Cauchy problem. We refer to [S] , [5] , [TU] for an overview 
of related results, which is complete at least for monotonous coefficients and provides sharp decay 
results for solutions. However, if the coefficient functions are allowed to bear a certain amount 
of oscillations, results may change dramatically. For the case of variable propagation speed this 
may even lead to exponentially growing energy, as pointed out in [12] or, at least destroy the 
structure of decay results, 0], [S], [S]. Recently, the first author developed a technique to obtain 
positive results for similar situations with strong oscillations by using a refined diagonalisation 
technique and a so-called stabilisation condition on coefficients, The aim of this note is to 
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extend this technique to the situation of oscillations in lower order terms, especially oscillations 
in dissipation terms. 

This note is organised as follows. At first we will introduce in Section [2] the assumptions we 
impose on the coefficient function b — b{t); particular examples of admissible coefficients are 
given in Section[31 The construction of £{t, s, ^) and will be done in Section^ where we introduce 
zones and give precise information on the structure of the fundamental solution depending on 
corresponding areas of the phase space. Finally, Section [5] collects the main results of this 
note. We present two theorems describing sharp energy decay results for solutions of the above 
introduced Cauchy problem. 

Throughout these notes we denote by c or C various constants which may change from line to 
line. Furthermore, f{p) < g{p) for two positive functions means that there exists a constant such 
that f{p) < Cg{p) for all values of the parameter p. Similarly, f{p) > g{p) means g{p) < f{p) 
and f{p) « g{p) stands for f{p) < g{p) and g{p) < f{p). For a matrix A we denote by ||A|| its 
spectral norm, while \A\ stands for the matrix composed of the absolute values of the entries. 



2. Assumptions 



Tools used in the approach are closely related to conditions on the coefficients. We impose that 
the dissipation can be written as 

where functions fi{t) and a{t) carry different kind of information: will determine the shape 
of the coefficient, while cr{t) contains oscillations (and zero mean in a certain sense). In detail 
our assumptions are 

(1) : fj.{t) > 0, ^'{t) < and \imsuptfi{t) < 1; 

t — 'OC 

(2) : generalised zero mean condition 



sup 

t 



a{s)ds 



(3): stabilisation condition 
exp 



a{s)ds — LOo 



< oo; 



< ei{t) = o{t) 



with a suitable (uniquely determined) lOoo > (and a function 0(t) normalised such that 
e(0) = /i(0)); 
(4): symbol-like conditions for derivatives 



< CfcS(t) 



1,2,.. 



with E{t) > Q{t); 
(5): together with the compatibility condition 



/>OC 



„\—rii— 1 



ds < e{ty 



between (3) and (4). 



Let's explain the philosophy behind the conditions (1) to (5): We assume that the reader is 
familiar with [9 and/or [2j. Since ^{t) should describe the shape of the coefficient, the assump- 
tions of (1) are related to The limit conditions excludes the exceptional case from [7], 
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where a structural change in the representation of solutions occurrs. Condition (2) describes 
that a(t) contains oscillations and the integral assumption implies a zero mean condition of cr(t). 
The stabilisation condition (3) is related to [2] and [3] (after a Liouville type transformation of 
variables). Note that (2) implies G(i) 0{t), the stabilisation improves this trivial estimate. 
We will use the notation 

/(t) a if and only if / |/(s) — a|ds = o{t) 

Jo 

for stabilising functions. Some elementary properties are given later on. The symbol- like esti- 
mates of assumption (4) are thought to be weaker than the ones from [8], [5, where — (l + t) 
was used. Stabilisation allows to use weaker assumptions on derivatives by shrinking the hyper- 
bolic zone to 

ZHyp{N) = {{t,o\em\>N}- 

We pay for this by using more steps of diagonalisation. The number of steps for diagonalisation 
will be the number m from condition (5). It implies that remainder terms are uniformly integrable 
over the hyperbolic zone after applying m steps. 

Remark 2.1. According to Appendix A. 6 and under assumption (2) the condition (3) holds if 
and only if (T(s)ds ^ logWoo. We will exploit this fact in the examples stated below. 



3. Examples 



We will collect some examples to illustrate the nature of our assumptions. 

Example 3.1. First we set /i(t) = with a fixed constant /i e (0, 1/2). Then (1) is fulfilled. 
Furthermore, (j{t) — /i(t)sin(i") satisfies (2) to (5). Indeed, (2) follows from 







* sin(s") , 1 
ds = — 

a 



1 + s 







0{a~l)/a ^ 



-de 



sin 9 



de'wSi(r) =0(1), 



while for (3) we use that the above integral converges by Leibniz criterium such that with 

1+s 



eMIo 



exp( 



ds) > 
sin(s") 



ds) 



< 



t 



to 

1-a 



1 + S 

sin(s") , 

ds 

1-f s 

= eit) = oit) 



d0 



de : 



exp 
sin(s) 



smis 



1 



ds 



dO < 



-ds) - 1 



-°'d9 



d0 



for a G (0,1). Furthermore, |dj^ sin(r)/(l + i)| < Ckil + Condition (4) and (5) are 
satisfied for m = 1 if we take E{t) = (1 + 1)^-°"^^. 

Example 3.2. Let /i(t) = j^. If we consider a{t) = /i(t) sin(i/ log(e + t)) we find (similarly) 

that (1) to (3) are satisfied with 8(t) = log^(e-l-i). Derivatives satisfy |d^ sin(t/ log(e + 1))| < 
Cfe(l -I- i)-i log(e + t)-^. We can not satisfy (4) and (5). 

Example 3.3. For (j{t) = (1 + t)^^ sa\[t°') assumptions (2) and (3) are satsified if q;,/3 > and 
1 < a + /3 < 2. We get Q{t) = t^-'^-fi and assumptions (4) and (5) can hold only if /3 > 1. 
Combining this with /i(t) = ^ < 1/2, we can choose m — 1 and S(t) — t(i+/3)/2-"/2_ 
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Example 3.4. If we consider ^{t) = {i+t)\lg{e+t) ' ^^^^ the function a{t) = ii{t) sm{t / log{e + t)) 
may be chosen. In this case &{t) = log(e + t). For (4) and (5) we choose to = 1 and S(t) = 
VT+llog(e + i), 

1 . 1 



s log^ s log i ' 

Example 3.5. We give a further example of a coefficient function with higher to. Let x G C'^(M) 
with suppx = [— Ijl]; lx(^)l < 1 and x(^)di = 0. For given sequences tj > 1 and 5j < 1, 
j — 1,2,..., of positive real numbers with tj + Sj^ < tj+i — ^J+i ctud a suitably chosen real 
number 7 > we define 



0, t^u^i, 



'3 "J • 

Then condition (2) holds because | /p cr(s)ds| < 2. For condition (3) we set Woo — 1 and estimate 
the integral as 

(T{s)ds dr < ^ 2 = 2n 

for f e [tn — S^i^, tn+i — S^+il- Thus with tn = n" assumption (3) holds for 2n = 2tl/°' « t^^"' — 
e{t) = o{t), i.e. if a > 1. We set J„ = t"'^ = Then n" + n"'' < (n + 1)" - (n + 1)"''' holds 

true for large n if 7 < 1. Derivatives satisfy ^l{t)a{t)\ < (1 + tjO-^C^+i) « {I + t)-'^^''+'^\ thus 
(4) holds with S(t) = (1 + i)^. Now we choose 7 such that (5) holds for a given number to, thus 

X} 

(l + r)-(™+i)^dT<t-"/", 

t 

i.e. (m + 1)7 > 1 and (m + 1)7— 1 = m/a. Hence we choose 7 = ^^^^^ + 7j^;jjqr[y, which is smaller 
than 1 for all to = 1,2,.... The function may be chosen as in Example 13. II or [3^ 

This example shows that for any given number to and any stabilisation rate Q{t) — t^/", a > 1, 
we find a coefficient 2b{t) = ^{t) + a{t) subject to (2)-(5). 

Remark 3.6. The results derived later on will show that essential influence on decay properties 
of solutions conies from the shape function /i(t), while the 'perturbation' a{t) has no influence 
on the decay rate at all. Note, that we do not require that a{t) is small corresponding to in 
some L°° sense. We only require, that the oscillations contained in a{t) are 'neatly arranged'. 



4. Construction of the fundamental solution £{t,s,£) 

Main point of our concern is how to use assumption (3) for small frequencies. We will introduce 
a bit more of notation and denote £^{t, s, ^) the fundamental solution for the case a{t) = 0, 

Properties of £fi{t, s,S) are studied in [S], [S] (as low-regularity theory, i.e. without using symbol 
classes and further steps of diagonalisation). The basic behaviour of 0, ^) can be summarised 
as follows: The phase space can be decomposed into three parts (two zones and subzones), 

• the dissipative zone z[f^^(N) ^ {{t,£) \ |^| < Nii{t)} and 

• the hyperbolic zone z\l^y^{N) — {(i,0 I ICI ^ divided into the regions where 
ICI > ^m(O) and where N ^i{t) < \^\ < iV^(O). 
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The latter subdivision is merely for convenience and does not stand for any deep structural differ- 
ences of the fundamental solution. The subdivision into zones is essential as the following results 
show. The constant N does not matter in this case. In the hyperbolic zone the fundamental 
solution behaves like £o(i, s,^) multiplied by A(s)/A(i) and a uniformly bounded and invertible 
matrix. The function 

1 



\{t) = exp J ^(s)ds 

contains (the essence of) the influence of dissipation. In contrast to that the fundamental solution 
0,^) behaves in the dissipative zone essentially like (see e.g. [T] for a neat argument) 



This bad behaviour (bad in the sense that it destroys the energy estimate appearing naturally 
within the hyperbolic zone) has to be compensated by assumptions on the data. Our assumption 
(1) on fi{t) implies 

||£^(t,0,Odiag(|C|/(0,l)ll < ^ 



with {0 = yrr^. 

4.1. Estimates in the dissipative zone. Following the argumentation of [T] or [5] we see that 
the fundamental solution £{t,0,^) satisfies within the dissipative zone 

Z,UN) = {it,0 I lei < N^^it)} = 4t(iV) 
the same estimates as sketched above for £p(t, 0,^). For later use we will denote the boundary 
of the dissipative zone by i^^''. 

Lemma 4.1. The fundamental solution £{t, 0, ^) satisfies for all (t, f) G Zdiss{N) the point-wise 
estimate 

i^«,o.oi<^(|^|:; \' 

Proof. The proof is mainly taken from [11, Lemma 2.1]. We rewrite the system as system of 
integral equations, denoting the entries of the rows of £{t,0,£^) as v{t,^) and w{t,^). This gives 

vit,0^m+^\C\ I «^(r,e)dr, 
Jo 

^(^'^) = T^^2 + i|ei-^ / A2(r)«(r,C)dr, 
A^(t) A'^(<) Jo 

where \{t) — exp(/p &(T)dr) and = (771,772) = (I7O) or 77 = (0,1) for the first and second 
column, respectively. 

We start by considering the first column. Plugging the second equation into the first one and 
interchanging the order of integration gives 

v{t, = 1- 1^1' t r >^^{0)v{9, Odddr 

Jo X-^[T) Jo 

Jo Je ^{t) 
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such that A^(t)|^|w(t, ^) satisfies an Volterra integral equation 

1^1 {t)v{t, = Ht, + f k{t, 9, {0H0, Ode 







with kernel k{t,6,^) = — I'fpA^Ci) /g dT/A^(T) and source term h{t,^) — |^|A^(t). Assumptions 
(1) and (2) imply h{t,£^) < 1 uniformly on Zdiss{N)- Representing the solution as Neumann 
series, 

Jo Jo 

implies the bound A^(t)|f ^ S L°°{Zdiss{N)) following from the kernel estimate 



(t,«)ez<i„,(A')io o<i<t</> (t.^je^d^.W Jo Je ^ C^) 



sup |fc(i,6i,0|d6' < sup |^pA^(4^^) 



« dr 



(1) 



= sup \i\'X\tf')f' T^dr<(|^|f)f <1, 

based on \k{t,0,^)\ « |CpA2(t) /g*dT/A2(r) (by (2)) and the monotonicity of t/X^it) for large 
time. The second integral equation implies the corresponding bound for w{t,^), 

m'it)\wit,0\<m f m'{r)\vir,0\dr<\^\ fdr<l 
Jo Jo 

uniformly on Zdiss{N). 

For the second column we use the same idea: Plugging the second equation into the first one 
yields the new source term |^| dT/A^(T) < 1/A^(t). Therefore the representation as Neumann 
series yields X'^{t)v{t,^) G L°°{Zdiss{N)) and integration with the second integral equation gives 
consequently X'^{t)w{t, E L°°{Zdiss{N)). The statement is proven. □ 

Corollary 4.2. The fundamental solution 0, S,) satisfies within the dissipative zone Zdiss{N) 
the norm- estimate 

||5(t,0,0diag(|e|/(0,l)|| <3^. 

Proof. It remains to estimate the entries in the first column of the product on the right. Lemma [4.1l 
gives the uniform bound (A2(t)|^||)-i|^|/(^) < 1/A2(t). □ 

Remark 4.1. Later on wc will only use the estimate 1/X{t) within this zone. 

4.2. Diagonalisation Estimates in the hyperbolic zone. We sketch the results, the idea 
of proof is essentially the same as in [2] or [5] . A Fourier multiplier a{t, ^ belongs to the symbol 
class Sf^{mi, 7712} if it satisfies the symbol estimate 

for all ^ ^hyp{N) and k < £ and all multi-indices a e Nq. These symbol classes satisfy 

natural calculus rules (like the ones from Proposition 6]). The condition E{t) > Q{t) gives 
embedding relations for these symbol classes of the form 

Sff{mi,m2} ^ Slj{mi + k,m2 - k}, k>0, 
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which is important for the application of the diagonahsation scheme sketched below. Assumption 
(4) guarantees that b{t) e iSJJrHOi 1}: such that for {t,£^) e Zhyp{N) the matrix A{t,£^) consists of 
a main part from 0} and the lower order corner entry from S^{0, 1} ^ S'^{1, 0} . 

The following diagonahsation scheme is merely standard and adapted from 2| . In an introductory 
step we diagonalise the homogeneous main part using AI = (J ^i), such that — M^^V 
satisfies for all (t,^) £ Zhyp{N) 

with Vo{t,^) = diagdCI +i6(t),-|^| + ib{t)) and i?o(i,C) = i&W(?o)- This system is diag- 
onal modulo Rq e S^{0, 1}. Now we apply an iterative procedure to diagonalise it modulo 
5™-'={-fc,fc+ 1}, k = l,...,m. 

Lemma 4.3. There exists a zone constant N > such that for any k — 0, 1, . . .m there exist 
matrices 

• Nk{t,£,) e 5]^^''{0,0}, invertible with inverse Nj^^ e 5]^'^''{0, 0} and tending to the 
identity as t oo for all fixed ^ ^ 0; 

. i?fc(t,0G5]^"'{-^,fc + i}; 

. Vkit,0 e 5™-'{l,0} diagonal, I?fc(t, ^ = diag(T+(i, 0, r^-^l^, 0) 
satisfying the operator identities 

(Dt -Vk- Rk)Nk = ^fc(Dt - Vk+i - Rk+i) 

for fc = 0, l,...,m— 1. 

Proof. The proof goes by direct construction. Assume for this that we have given a system 
Dtl^W = withl?fc(t,0 = diag(T+(t,e),r-(t,a) € 5™^'{1,0} satisfying 

l4(t,0l-k,+(i,e)-T,7(t,0l>t^fclCI 

and antidiagonal Rk{t,S,) G S^^^''{ — k,k + 1}. Then we denote the difference of the diagonal 
entries as Sk{t,£^) = Tj^(t,^) — T^^{t,£,) and set 

-{Rk)i2/Sk\ 
{Rk)2i/5k ) 

such that [Dk^Nk] — —Rk and therefore 

= (D, _ - Rt,)Nk - Nk{Bt - Vk) 

= DtNk-[Dk,Nk]-RkNk 

= [BtNk) - RkiNk - /) e 5™-'=-i'°°{-fc -l,k + 2}. 

The matrix Nk{t,£,) is invertible, if we choose the zone constant N sufficiently large. This 
follows from the symbol estimate / — Nk € S^^''{—k — 1, fc + 1}. Thus by defining V^+i = 
Vk - diag(iVj:^B('=+i)) and Rk+i = diag(iVjr^B('^+i)) - N^^B^^+^'^ we obtain the operator 
equation 

(Dt -Vk- Rk)Nk = NkiBt - Vk+i - Rk+i) 
and it is easily checked that the assumptions we made are satisfied again. □ 

Finally we obtain for k = m that the remainder Rm{t, ^ 'S%{—'m, m+1} is uniformly integrable 
over the hyperbolic zone, 

/•OO /"OO 

/ \\RUt,0\\dt < i^r™ / sw-™-Mt < i^r^eif V™ < N, 



Nkit,0=I + 
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where we defined as the maximum of and the imphcitly defined zone 

boundary e(t^ OICI = 

N. To complete the construction of our representation we need more information on the diagonal 
matrices T>k- 

Lemma 4.4. For all k ^ 0,1, ... ,m the difference of the diagonal entries of Vkit, ^) is real. 

Proof. Again we proceed by induction over k and follow the diagonalisation scheme. For fc = 
the assertion is satisfied and the hypothesis 

(H) Rk{t, S,) has the form Rk = i( ^ ^'') with complex- valued f3k{t, £,) 

is true. Thus the construction imphes Nk = ^ + ^( ) ^itl^ det Nk — I — \Pk\'^/Sl ^ (for 

our choice of the zone constant N). Following [2] we obtain (with = |/3fcP/(5fc) 

N^^iVk + Rk)Nk = Y^^( diag (t+ - dkT+ - 5kdk,T^ - dkT^ + Skdk) + dkRk) 

and 



1 Uf -^^1 



such that Re^dt^ = ^ = Re^dt§^ imphes 



<5fc <5fc - 2 - -^"-^ 5fc 

,± _ ±_ 1 f J . T.^ f l^k ^ PkW : (^tdk 



Hence dk+i is real again and Rk+i satisfies (H) and, therefore, both statements are true for all 
k up to m. □ 

Now the construction of the fundamental solution £{t, s, ^) is merely standard. At first we solve 
the diagonal system Dt — I?m(t,^). Its fundamental solution is given by 



exp(^i^ V^i9,Od9^ =diag(e'/>' 



+ (e.«)de^gi/V;;(e.?)de 



Since 6m = — is real, it follows that Imr+ = Imr„j ImT™ and thus the matrix 



exp(^^ lmTmiO,Odd^ exp (^i Vm{e,(,)d6 



is unitary. Note, that the first factor is scalar. Now the integrability of the remainder term 
Rm{t, £,) over the hyperbolic zone implies (like in [1]) that the fundamental matrix of Dt — Vm — 
Rm is given by exp(i 'Dm{6,£_)d6)Q„i{t,s,£,) with a uniformly bounded and invertible matrix 

Qm.{t,S,£,), 



Qmit,S,0^I+J2 I Rm{h,S,C) I Rm{t2,S,i)--- j (tfc , S, Odtfc ' ' ' dtsdii 

where 

R,nit,s,0=exp{-i f Vmi0,Ode)Rm{t,OeM^ [ 'Dra{6,i)de) 

J s J s 

is an auxiliary function. The representation implies that Qm(^, tends to the identity as 
i — > oo locally uniform in s and ^. Collecting these results we obtain 
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Lemma 4.5. The fundamental matrix £{t,s,^) satisfies 

£it,s,0=M-' (j[ N-\t,0^ exp (^iJ^Vm{e,^)de^ Qm{t,s,^) ^ Nk{s,0^ M 

for all {t,£,), (s,^) G Zhyp{N), where 

• the matrices Nk(t,£,) are uniformly bounded and invertible with Nk{t,£,) —^ I and 

• Qmit,s,^) is uniformly bounded satisfying Qm(t,s,£,) Qm{oo,s,S,), 

both limits locally uniform in ^ as t ^ oo. 

Hence, the time-asymptotics of solutions is encoded in the real- valued function Imrm(i,^), 



||£:(t,s,OII «exp(^-^*ImT™(^,Od^) , t 



oo 



locally uniform in ^ for fixed s (and such that (s, ^) G Zhyp{N)). We can use our representation 
of Tm{t,^) to deduce 

m— 1 Q , 

such that 



— ' A(£) 
~ x{t)- 



Corollary 4.6. The fundamental matrix £{t,s,^) satisfies uniformly in (s,^) G Zhyp{N) 
the two-sided estimate 

4.3. Estimates in the intermediate zone. Since Q{t)n{t) < 8(/:)/(l+t) = o(l), there remains 
a gap between the dissipative and the hyperbolic zone. We will denote this zone as 

Zint{N) = {{t,^)\tf^ <t <tf}. 

In this zone we relate £{t,s,£) to (known estimates for) £n{t,s,£,) and use the stabilisation 
condition (3). For this we solve 

DtA{t, s, = {^{t, - Mt^ 0)Ht, s, 0, A(s, s, = I, 

which gives 

A(t,5,c) = diag(^i,cxp(-^ <7{e)de)y 

and make the ansatz £{t, s, ^) = A{t, s, £,)TZ{t, s, ^). It follows that the matrix TZ{t, s, ^) satisfies 

Dtn{t,s,o ^ A{s,t,OAf,{t,OHt,s,on{t,s,o, n{s,s,o = i, 

where the coefficient matrix in this system has the form 
A^{t, s, = Ms, t, OA^{t, OA{t, s, 

exp(-/>(0)d0)|^|\ 
l^exp(/>(0)d0)|C| Ht) ) 

Note that condition (3) means A{t, s,^) diag(l, ti)oo(s)^^), where we use tdoo{s) = Wqo exp(— a{6)d9). 
Condition (2) implies that < c < a>oo(s) < C < oo with suitable constants. Thus the new 
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speed of propagation satisfies the stabilisation condition (as used in l2j), while the dissipation 
term has no bad influence as consequence of assumption (f). 

We denote by Afj,{t, s, ^) the matrix 



iOoo{s)-^\^\ 



and solve the corresponding system Dt — A^^. The diagonaliser of the |^|-homogeneous part is 
given by 

f -1 



M{s) = 







1 1 



such that 
Surprisingly, this means 

MM-\s)A^[t,s,OM{s)M-^ = A^{t,0, 
such that the solution s,^) to the auxiliary problem 

(Dt - 4(t, S, - 0, S, - /, 

satisfies s,^) = M(s)Af s, ^)MM-i(s). This relation implies: 

Lemma 4.7. The matrix £^{t, s,C) satisfies uniformly in G ^hypi^) 



A(£) 

m 



Now we use the stabilisation property of s,^) to find TZ{t,s,S,) of the form TZ{t,s,£,) 

^)Q7j(i, s, ^). The coefficient matrix of the differential equation satisfied by Qn, 

^tQnit, s, = t, OiMt, - OQK(t, s, 0, Qk(s, s, = /, 

satisfies the estimate (note, that the two-sided estimates for £^ cancel each other) 

||4(s, r, - 4(r, 0)4 (r, Olldr 



exp 



tT(r)dr - Woo(s) 



del 



y exp ^ tT(T)dr^ - 



d6l < |C|e(t) < N. 



Now the standard construction of QK(i, s,^ in terms of the Peano-Baker series (as done for 
Qm{t, s,^) before) gives uniform bounds for this matrix and for its inverse within the intermediate 
zone. Thus 

Lemma 4.8. The fundamental matrix £{t, s,S^) can he represented in Zint{N) as 

£{t, s, i) = A{t, s, 04 (i, 5, OQnit, s, 0, 

where A{t,s,^) and Qiz{t,s,^) are uniformly bounded in {t,£,),{s,^) G Zi„t(7V) and 4(^1^,0 
satisfies the bound of Lemma \4. 7[ 
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Corollary 4.9. The fundamental matrix £{t, s,^) satisfies uniformly in {t,^),{s,^) € Zint{N) 
the two-sided estimate 

5. Main results - Energy decay estimates 

The results of the previous section can be collected as energy estimates for solutions to the 
original Cauchy problem. We obtain two results. The first one states the energy decay estimate. 
The use of the iJ^-norm for the first datum is essential for the validity of the statement. This 
cancels the |^|^^ in the estimate of Lemma [4. II 

Theorem 1. Solutions to the Cauchy problem 

au + 2b{t)ut^0, U(0, Ut(0, •)=M2 

for data ui e H^{M") and U2 G L^(R") and coefficient function b(t) subject to conditions (l)-(5) 
satisfy the a-priori estimate 

||Vu(t, + \\ut{t, < C^{\\ui\\hi + \\u2\\l-) 

with a constant C depending on the size of the coefficient b{t) and its first m derivatives, where 
the function \{t) is given by \{t) = exp(l/2 /i(T)dr) « exp(J'p 6(s)ds. 

Proof. By Plancherel's theorem it is equivalent to prove the corresponding statement in Fourier 
space, 

+ \\ut{t, < C^{\mui\\L- + \\u2\\l-). 

which reduces by the unitarity of Riesz transform, ||^u||i2 = || |(^|u||i2, and in view of our system 
reformulation to the estimate 

||£(t,0,Odiag(|e|/(e>,l)K,||<^||K,|| 

for all Vb e and uniform in ^. But this is just the combination of Corollaries 14.21 14.61 
and [131 □ 



The second result is an application of Banach-Steinhaus theorem on the (dense) subspace of data 
for which does not belong to the Fourier support. It follows essentially from the fact that the 
matrix s, ^) tends locally uniform to a invertible matrix Q(cxd, s, ^) inside the hyperbolic zone 
Zhyp{N). The proof is analogous to the corresponding one in 9, Theorem 31], TT, Corollary 
3.2]. 

Theorem 2. For any fixed choice of data ui G H^{W") and U2 G -L^(R") we find constants c and 
C such that under the assumptions of Theorem 1 the solution to the Cauchy problem satisfies 

c<\{t){\\Vu{t,-)\\L2 + \\ut{t,-)\\L2)<C. 

Proof. The proof consists of two parts. 
Part 1. We denote 

'M-^\-Ht^) exp (i //^ Vra{e,£.)de) M, t > t^, 
x-^{t)i, t<t^, 
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(2) 

where we used for convenience the abbreviation t^^ ~ . In a first step we show that the limit 

t — >oo 

exists uniformly on |^| > c for any given c > and defines an invertible matrix W{^) for 7^ 0. 
Indeed, the representation of Lemma 14.51 shows that 



M-^exp l^-i^ Vmie,Ode^ M£it,^,0 Af-iQ„,(oo,t^,0 (^J[ Nk^^^^ M, 

where the difference between left and right hand side can be estimated by ((1 which 
tends uniformly to zero on any set |^| > c with c > 0. Note, that the right hand side is 
independent of c and belongs to L°°(R") and W{^) is obtained after multiplication by £{t^,0, ^). 

Therefore, W{£,) € L°°(R") is well-defined and it remains to check the invertibility. For this we 
apply Liouvillc theorem to our initial system. This gives 



det£(i,s,^) ^ cxp ^i^ trA{T,()dT^ = exp (^^2^ fe(T)dr 



while 

detexp l^-i^ Vmie,Ode^ = exp (2 T„(0,C)d0^ « 

Therefore, | detyV(^)| w 1 and thus W{^) and W^^{^) are both uniformly bounded. 

Note, that the matrix ^) is a scalar multiple of a unitary matrix, the factor essentially given 
byA-i(t). 

Part 2. We consider the dense subspace 

L|5I>, = {./ e : dist(0,supp/) > c} C L^{R'';C'^). 

In Part 1 we have shown that £^^{t,^)£{t,0,^) W(^) uniformly in |^| > c. Therefore, on the 
operator-level 

C^(i,D)£(t,0,D) W(D) 
pointwise on L^^|^^. Theorem 1 provides us with the norm-estimate 

||f(i,0,D)diag(|D|/(D),l)|U2^i2 <1, 
such that by Banach-Steinhaus theorem the strong convergence 

s-\im£-\t, B)£{t, 0, D) diag(|D|/(D), 1) = >V(D) diag(|D|/(D), 1) 

t — *oo 

follows on L2(R";C2). Therefore, for all Vq G diag(|D|/(D), 1)L^{W\C^) 

\\£-\t,B)£{t,0,B)Vo ^Wi-D)Vo\\L- 
^\\\{t)£it,0,B)Vo - m£4t,D)W{D)Vo\\L- -> 0, 

while we already know that ||A(i)£*(i,D)W(D)Vb||L2 w 1. Therefore, the first term satisfies 
two-sided bounds 0, D)Vb||L2 sa ^^'^{t) and the theorem is proven. □ 
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Appendix A. Elementary properties of stabilising functions 

We collect some properties of stabilising functions. For / e Lj^^{R+) and a e M we write f a 
iijl\f{s)-a\ds = o{t). 

1. a is uniquely determined. Assume the condition is also satisfied for q', then 

t\a-a'\= I \a-a'\ds< [ \f{s)-a\ds+ [ \f{s)-a'\ds = o{t) 
Jo Jo Jo 

and the assertion follows. 

2. The number a does not depend on the lower bound of the integral, the number is used for 

convenience only. 

3. If /i ^ ai and /2 ^ a2 then /i + c/2 ^ ai + ca2. 

4. Stabilisation f a does not imply convergence of f{t), but if we additionally know that the 
limit exists, then it must be equal to a. 

5. If / ^ a and g is monotone with g'{t) bounded, then f o g a. This follows by substitution 
in the integral, 

r9'\t) ft 

/ \f{g{s'))-a\ds'= \f{s)-a\ds = o{t)=o{g-\t)). 

Jg-HO) Jo 

6. On the contrary, if / ^ q and g is Lipschitz, then g o f g{ct)- Indeed, the Lipschitz 
condition implies directly 

* Wis)) - 9{a)\ds <L f \f{s) - a\ds = o{t). 
Jo 

Note that it is sufficient to require that g is locally Lipschitz if / is bounded, i.e. we need the 
Lipschitz condition on the range of /. 

Acknowledgements. These notes are based on discussions on the subject held when the second 
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foundation (DFG) as project 446 JAP 111/3/06. 



References 

[1] J. Emmcrling. Wave equations with timc-dcpcndcnt coefficients. Adv. Math. Sci. Appl, to appear. 

[2] F. Hirosawa. On the asymptotic behavior of the energy for wave equations with tinie-depcnding coefficients. 

Math. Arm., 3.39(4) :819-839, 2007. 
[3] F. Hirosawa, J. Wirth. Generahzed energy conservation law for the wave equations with variable propagation 

speed, in preparation. 

[4] M. Reissig and K. Yagdjian. One application of Floquet's theory to Lp-Lq estimates for hyperbolic equations 

with very fast oscillations. Math. Methods Appl. Sci., 22(11);937-951, 1999. 
[5] M. Reissig and K. Yagdjian. About the influence of oscillations on Strichartz-type decay estimates. Rend. 

Semin. Mat., Torino, 58(3):375-388, 2000. 
[6] M. Reissig and J. Smith. L^-L^ estimate for wave equation with bounded time dependent coefficient. Hokkaido 

Math. J., 34(3):541-586, 2005. 
[7] J. Wirth. Solution representations for a wave equation with weak dissipation. Math. Meth. Appl. Sci., 27:101- 

124, 2004. 

[8] J. Wirth. Asymptotic properties of solutions to wave equations with time- dependent dissipation. PhD thesis, 
TU Bergakademie Freiberg, 2005. 

[9] J. Wirth. Wave equations with time-dependent dissipation I. Non-effective dissipation. J. Differential Equa- 
tions, 222:487 514, 2006. 

[10] J. Wirth. Wave equations with time-dependent dissipation. II: Effective dissipation. J. Differential Equations, 
232(1):74-103, 2007. 



14 



FUMIHIKO HIROSAWA AND JENS WIRTH 



[11] J. Wirth. Scattering and modified scattering for abstract wave equations with time-dependent dissipation. 

Adv. Differential Equations, 12(10):1115-1133, 2007. 
[12] K. Yagdjian. Paramctcric resonance and nonexistence of the global solution to nonlinear wave equations. J. 

Math. Anal. Appi, 260(l);251-268, 2001. 

F. HiROSAWA. Department of Mathematical Sciences, Yamaguchi University, 753-8512, JAPAN 
E-mail address: hirosawa@yamaguchi-u.ac.jp 

J. Wirth, Department of Mathematics, Imperial College, London SW7 2AZ, UK 
E-mail address: j.wirth@imperial.ac.uk 



